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Transition State Theory overestimates reaction rates in solution because conventional dividing 
surfaces between reagents and products are crossed many times by the same reactive trajectory. 
We describe a recipe for constructing a time-dependent dividing surface free of such recrossings in 
the presence of noise. The no-recrossing limit of Transition State Theory thus becomes generally 
available for the description of reactions in a fluctuating environment. 
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Chemical reaction rates are determined by what hap- 
pens in a very small fraction of phase space: Such is 
the premise of Transition State Theory (TST). It re- 
volves around the concept of the Transition State (TS), 
also known as the activated complex, that is defined by 
the dividing surface that best separates "reactant" and 
"product" regions in phase space. Actually, the TS is 
a general structure found in all dynamical systems that 
evolve from "reactants" to ''products". In addition to 
chemical reaction dynamics ^ also determines rates 
in many other interesting systems such as asteroid cap- 
ture , mass transport through the solar system , the 
rearrangements of clusters Q , the ionization of atoms || , 
conductance due to ballistic electron transport through 
microjunctions Q, and diffusion jumps in solids Q. 

TST is central to the understanding of chemical re- 
actions because its radical simplification captures the 
physics involved in a pictorial way, but its accuracy 
hinges on the so-called "no-recrossing" requirement: If 
reagents and products are separated by an abstract hy- 
persurface in phase space, the rate will be exact if each 
reacting trajectory crosses this surface once and only 
once. Recrossings will lead to an overestimate of the 
rate, aproblem that occupied Wigner in the early days of 
TST la, |2j ■ While in low-dimensional systems a solution 
has been known for a long time [lfl ]. enforcing the no- 
recrossing condition in more than two degrees of freedom 
has been an open problem until very recently [ill El fl3| . 
It is even more difficult to uphold for reactions in ran- 
domly fluctuating environments (such as chemical and 
biological processes in liquids) where the noise causes the 
particle to move randomly back and forth, so that it will 
typically cross any fixed dividing surface in phase space 
many times Eil- To solve this long-standing problem, 
many approximate TS structures have been suggested 
in the literature El El El El El El E3 Martens 
|2lj |. in particular, constructed a stationary stochastic 
separatrix (the collection of all phase-space points for 
which the reaction probability is 1/2) that foreshadows 
the time-dependent structures to be described here. In 
some cases, infinite dimensional representations of the 
Langevin equation El El EH; in which the system is 
coupled to a bath of harmonic oscillators, have been used 



to obtain approximate infinite-dimensional dividing sur- 
faces that lead to excellent approximations to the rate 

E1E1E3- 

In this Letter, an exact transition state in a noisy envi- 
ronment is constructed as a dividing surface which moves 
so as to avoid recrossings. For each instance of the noise, 
there is a unique trajectory which remains in the vicin- 
ity of the barrier for all time. This trajectory serves as 
the origin of a moving coordinate system to which the 
geometric structures of the noiseless phase space are at- 
tached. The result is a moving dividing surface, free of 
recrossings, that describes a reaction influenced by noise 
in the same way as a static dividing surface does in con- 
ventional TST. This construction goes beyond the geo- 
metric approach of E3| by fully taking into account the 
time-dependence of the fluctuating force. Indeed, the 
moving separatrices introduced here determine with cer- 
tainty if a trajectory will or will not react. 

We illustrate our findings on a reactive system de- 
scribed by the Langevin equation Ef| 



(1) 



The vector q denotes a set of N mass-weighted coordi- 
nates, U(q) the potential of mean force governing the 
reaction, T a symmetric positive-definite friction matrix 
and £ a (t) a fluctuating force assumed to be Gaussian 
with zero mean. The subscript a represents randomness 
by labeling different instances of the fluctuating force. 
The latter is related to the friction matrix T by the 
fluctuation-dissipation theorem |2f| 



£«(*)&(0) =2k B Trs(t-t') 
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where the angular brackets denote the average over the 
instances a of the noise. 

The reactant and product regions in configuration 
space are separated by a potential barrier whose position 
is marked by a saddle point q^ = of the potential U (q) . 
The reaction rate is primarily determined by the dynam- 
ics in a small neighborhood of El El El El El El , 
so that the deterministic force can be linearized around 
the saddle point to obtain 



t a (t) = nq a (t) - TqT a {t) + &(t) 



(3) 
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with a symmetric matrix Jl given by fijj = 
— (d 2 U I 'dqidqj) -. -4 ■ Coordinates can always be chosen 
to that fl is diagonal. It must have one positive eigen- 
value equal to the squared barrier frequency and N — 1 
negative eigenvalues — Lof, where u>i are the frequencies 
of transverse oscillations. Including velocity-dependent 
forces from, e.g., a magnetic field [5| in this setup would 
be straightforward. Note that no restriction on the di- 
mensionality of the configuration space is imposed. 

Our aim in the following is to construct a recrossing- 
free surface for equation J2J. The difference 

Aflt) = £(i)-g*(t) (4) 

of two trajectories under the influence of the same fluc- 
tuating force, which describes the location of the trajec- 
tory q a (i) relative to the moving origin q%(t), satisfies a 
deterministic equation of motion. Since this difference 
is independent of noise, there is no need to refer to a 
in Aq(t). Thus, invariant manifolds and a no-recrossing 
surface can be specified for the relative dynamics. These 
geometric objects in the noiseless phase space can then 
be regarded as being attached to the reference trajec- 
tory q%(t) while being carried around by it. They define 
randomly moving invariant manifolds and a moving no- 
recrossing surface in the phase space of the original, noisy 
system. 

The construction outlined above can be carried out for 
any reference trajectory q\{t) and leads to a multitude 
of no-recrossing surfaces, but only one specific surface is 
relevant to the reaction dynamics: The crossing of this 
surface should signal the transition of the trajectory q a (t) 
from the reactant side to the product side of the bar- 
rier. If the reference trajectory is chosen arbitrarily, it 
will typically descend into either the reactant or prod- 
uct wells over time. However, only a reference trajec- 
tory that remains in the vicinity of the barrier can carry 
a no-recrossing surface that actually describes the reac- 
tion. Indeed, as shown here, for each instance of the noise 
there is a unique reference trajectory with this property. 
It represents, in mathematical terms, an invariant mea- 
sure of the noisy dynamical system 27] . We call it the 
Transition State Trajectory and henceforth restrict the 
notation q\{t) to this privileged reference trajectory. 

We solve equation Q by rewriting it as a first-order 
equation in 2iV-dimensional phase space with coordinates 
x = (q,v), with v — q, and diagonalizing its linear part. 
Equation then decomposes into a set of 2N indepen- 
dent scalar equations 

%aj{t) — £j x aj{t) + £,aj(t) j (5) 

where ej are the eigenvalues of the linear part, x a j the 
components of x in a basis of eigenvectors and £ Q j the 
corresponding components of (0,£ Q (t)). 

The components Aqj(t) of the relative coordinate (0J 
satisfy a noiseless version of (0 with £ a j = 0. They grow 
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FIG. 1: Phase portrait of the one-dimensional noiseless 
damped dynamics corresponding to equation @. The thick 
line indicates a possible choice for the no-recrossing surface. 
The simple relative dynamics illustrated here is exemplified 
by the trajectories in figure Ufa). 

or decay exponentially, depending on whether the eigen- 
value ej has positive or negative real part. In N = 1 de- 
gree of freedom there is one positive and one negative real 
eigenvalue. The phase portrait of the dynamics is shown 
in figure ^ The eigenvectors span one-dimensional sta- 
ble and unstable manifolds of the saddle point. They act 
as separatrices between reactive and non-reactive trajec- 
tories. The knowledge of the invariant manifolds allows 
one to determine the ultimate fate of a specific trajec- 
tory from its initial conditions. It is also easy to identify 
lines in the quadrant of reactive trajectories that are sur- 
faces of no recrossing. Clearly, the half- line Aq = 0, with 
Aq > serves this purpose. 

In multiple degrees of freedom, transverse damped os- 
cillations must be added to the phase portrait in figure^ 
Their presence manifests itself through N — 1 complex 
conjugate pairs of eigenvalues e,-. For strong damping, 
some of the transverse modes can become overdamped, 
so that further eigenvalues become negative real. In any 
case, there is exactly one positive real eigenvalue which 
corresponds to the particle sliding down the barrier. In 
all other directions, the dynamics is stable. The eigen- 
vector corresponding to the smallest negative eigenvalue 
together with the unstable eigenvector span a plane in 
phase space in which the dynamics is given by the phase 
portrait of figure ^ The separatrices between reactive 
and nonreactive trajectories that were identified for the 
one-dimensional dynamics together with the stable trans- 
verse subspace form separatrices in the high-dimensional 
phase space. In a similar manner, a no-recrossing sur- 
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FIG. 2: A random instance of the TS trajectory (black) and a reactive trajectory (red) under the influence of the same noise 
in a system with N = 2 degrees of freedom: the unstable reactive coordinate q u and the stable transverse coordinate q s - The 
potential is U(q u , q B ) = — k^bQu + f^s^s • Trajectories are projected onto (a) the reactive degree of freedom, (b) the transverse 
degree of freedom, and (c) configuration space. Units are chosen so that Wb = 1 and feT = 1. The transverse frequency is 
uj s = 1.5, and the friction is isotropic, T — 7I, with 7 = 0.2. The bottom of each column shows the projected trajectories in 
the corresponding space. Above this, their time evolution is illustrated using the same axes. The blue cut marks the unique 
reaction time £ rea ct = 8.936, when the moving TS surface is crossed. Dotted lines in this cut, at t = 0, and at t — 20 indicate 
the moving coordinate axes centered on the TS Trajectory. These axes are labeled explicitly only at the top face. Dashed lines 
in the cuts of column (a) show the moving invariant manifolds. No TS is indicated in column (b) because the q a — v s subspace 
lies entirely within the moving TS surface. Thick green dots indicate repeated crossings of the stationary TS surface q u — 0. 
Green lines in the top face of all three columns show the reactive trajectory in relative coordinates Aq, Av (not to scale in (a), 
for graphical reasons). The typical behavior of a reactive trajectory shown in figure is visible in the top face of column (a). 



face in the full phase space is spanned by a no-recrossing 
curve in the plane and the transverse directions. 

Using the Green function technique and assuming 
3tej < 0, we find a particular solution 

<*•(*) = / e-^£ aj (t + T)dT (6) 
J — 00 

of equation JSJ . The general solution consists of © plus 
an exponential term that must be suppressed to keep the 
solution bounded as t — -> —00. Similarly, for unstable 
components with dtej > 0, the only solution to J5j that 
remains bounded for t — > +00 is 

/>oo 

xl){t) = - e-^^(t + r)dr. (7) 
Jo 

Equations © and (7J specify all the phase space com- 



ponents of the unique random trajectory — viz. the TS 
Trajectory — that remain bounded in the remote future 
as well as in the distant past. Its configuration space rep- 
resentation is obtained by transforming back to position 
and velocity coordinates. Because the components £ QJ - of 
the fluctuating force are Gaussian random variables with 
zero mean, so are the components (0 and J7J) of the TS 
trajectory. Their joint probability distribution can read- 
ily be specified through the explicit calculation of their 
cross-correlation function by insertion of the correlation 
function of the noise. 

As described above, the TS trajectory serves as the ori- 
gin of a moving coordinate system to which the geometric 
structures in the noiseless phase space are attached. This 
is illustrated in figure |21 for = 2 degrees of freedom. 
The phase space is four-dimensional and already exhibits 



4 



all the salient features of the dynamics in arbitrarily high 
dimensions. The figure compares a TS trajectory with a 
reactive trajectory under the influence of the same noise. 
The trajectory approaches the TS trajectory from the re- 
actant side q u < 0, remains in its vicinity for a while and 
then wanders off to the product side. Because the mov- 
ing invariant manifolds of the TS trajectory are known, 
it can be predicted already at time t — that the trajec- 
tory actually will lead to a reaction instead of returning 
to the reagent side of the saddle. 

The reactive trajectory crosses the space-fixed TS sur- 
face q u = several times before it finally descends on 
the product side. The moving TS surface Aq u = 0, by 
contrast, is crossed only once, at the reaction time t rcact 
indicated by blue lines. That this is actually the case can 
be seen from the curves in the top faces of each column, 
which indicate the noiseless relative motion Aq(t), Av(t) 
between the TS trajectory and the reactive trajectory. 
The relative coordinates show, as expected, the hyper- 
bolic motion known from figure ^hi the unstable degree 
of freedom, a damped oscillation in the stable transverse 
degree of freedom, and a superposition of the two in the 
configuration-space projection in figure Ofc). 

The derivations presented here are based on the as- 
sumption (J2J of white noise. If the noise were corre- 
lated, the friction would contain memory and the dy- 
namics would consequently be significantly more compli- 
cated 26] . The construction of the TS trajectory can 
be carried out nonetheless, ft leads to integral formulas 
similar to © and 0, with the eigenvalues ej given by 
the Grote-Hynes equation 18, 28]. These generalizations 
will be presented in a forthcoming publication. 

In recent simulations [2j| [3(| , ensembles of transition 
paths, i.e., trajectories connecting reactants to products, 
have been shown to be useful in computing rates even 
for high-dimensional systems. Our identification of the 
TS trajectory provides additional insight into the geome- 
try of reaction dynamics that complements the transition 
path ensemble method. While the TS trajectory intro- 
duced here is not itself a transition path, it carries a 
dividing surface that is transverse to all transition paths 
and is crossed once and only once by each of them. It is 
no harder to compute than a typical trajectory. This is to 
be contrasted with most TST approaches in which find- 
ing a no-recrossing dividing surface requires the onerous 
solution of a highly nonlinear problem. 

In summary, to generalize the formalism of TST to 
reactive systems driven by noise, a moving surface has 
been constructed that is crossed once and only once by 
each transition path. This surface is therefore suited to 
take over the role of the well-known and widely applied 
static dividing surface of conventional TST in a broader 
time-dependent setting. 



This work was partly supported by the US National 
Science Foundation and by the Alexander von Humboldt- 
Foundation. 



[1] 
[2] 

[3] 

[4] 

[5] 

[6] 
[7] 

[8 
[9 
[10 

[11 

[12 

[13 

[14 
[15 

[16 

[17 
[is; 

[19 

[20 

[21 
[22 
[23 

[24 

[25 
[26 

[27 

[28 

[29 

[30 



W. H. Miller, Faraday Discuss. Chem. Soc. 110, 1 (1998). 
C. Jaffe, S. D. Ross, M. W. Lo, J. E. Marsden, D. Farrelly, 
and T. Uzer, Phys. Rev. Lett. 89, 11101 (2002). 
W. S. Koon, M. W. Lo, J. E. Marsden, and S. D. Ross, 
Chaos 10, 427 (2000). 

T. Komatsuzaki and R. S. Berry, Adv. Chem. Phys. 123, 
79 (2002). 

C. Jaffe, D. Farrelly, and T. Uzer, Phys. Rev. Lett. 84, 
610 (2000); Phys. Rev. A 60, 3833 (1999). 

B. Eckhardt, J. Phys. A 28, 3469 (1995). 

M. Toller, G. Jacucci, G. DeLorenzi, and C. P. Flynn, 
Phys. Rev. B 32, 2082 (1985). 

E. P. Wigner, Faraday Discuss. Chem. Soc. 34, 29 (1938). 

E. P. Wigner, J. Chem. Phys. 7, 646 (1939). 

P. Pechukas and F. J. Lafferty, J. Chem. Phys. 58, 1622 

(1973). 

T. Komatsuzaki and R. S. Berry, J. Chem. Phys. 110, 
9160 (1999). 

S. Wiggins, L. Wiesenfeld, C. Jaffe, and T. Uzer, Phys. 
Rev. Lett. 86, 5478 (2001). 

T. Uzer, C. Jaffe, J. Palacian, P. Yanguas, and S. Wig- 
gins, Nonlinearity 15, 957 (2002). 

E. Pollak and P. Talkner, Phys. Rev. E 51, 1868 (1995). 
J. C. Keck, Adv. Chem. Phys. 13, 85 (1967). 
J. T. Hynes, Theory of Chemical Reaction Dynamics, 
vol. 4 (CRC, Boca Raton, FL, 1985). 

D. G. Truhlar, B. C. Garrett, and S. J. Klippenstein, J. 
Phys. Chem. 100, 12771 (1996). 

E. Pollak, J. Chem. Phys. 85, 865 (1986). 

S. C. Tucker, in New Trends in Kramers' Reaction Rate 
Theory, edited by P. Hanggi and P. Talkner (Kluwer Aca- 
demic, The Netherlands, 1995), pp. 5-46. 
E. Pollak, in Theoretical Methods in Condensed Phase 
Chemistry, edited by S. D. Schwartz (Kluwer Academic 
Publishers, Dordrecht, 2000), pp. 1-46. 

C. C. Martens, J. Chem. Phys. 116, 2516 (2002). 
R. Zwanzig, J. Stat. Phys. 9, 215 (1973). 

A. O. Caldeira and A. J. Leggett, Phys. Rev. Lett. 46, 
211 (1981), Ann. Phys. (N. Y.) 149, 374 (1983). 
E. Pollak, H. Grabert, and P. Hanggi, J. Chem. Phys. 
91, 4073 (1989). 

R. Graham, J. Stat. Phys. 60, 675 (1992). 
R. Zwanzig, Noneqmlibrium Statistical Mechanics (Ox- 
ford University Press, London, 2001). 
L. Arnold, Random Dynamical Systems (Springer, 
Berlin, 1998). 

R. F. Grote and J. T. Hynes, J. Chem. Phys. 73, 2715 
(1980). 

C. Dellago, P. Bolhuis, F. S. Csajka, and D. Chandler, J. 
Chem. Phys. 108, 1964 (1998). 

P. G. Bolhuis, D. Chandler, C. Dellago, and P. Geissler, 
Annu. Rev. Phys. Chem. 53, 291 (2002). 



